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Abstract 

Given a family of feasible subsets of a ground set, the packing prob¬ 
lem is to find a largest subfamily of pairwise disjoint family members. 
Non-approximability renders heuristics attractive viable options, while 
efficient methods with worst-case guarantee are a key concern in computa¬ 
tional complexity. This work proposes a novel near-Boolean optimization 
method relying on a polynomial multilinear form with variables ranging 
each in a high-dimensional unit simplex, rather than in the unit interval 
as usual. The variables are the elements of the ground set, and distribute 
each a unit membership over those feasible subsets where they are in¬ 
cluded. The given problem is thus translated into a continuous version 
where the objective is to maximize a function taking values on collections 
of points in a unit hypercube. Maximizers are shown to always include 
collections of hypercube disjoint vertices, i.e. partitions of the ground 
set, which constitute feasible solutions for the original discrete version of 
the problem. A gradient-based local search in the expanded continuous 
domain is designed. Approximations with polynomial multilinear form of 
bounded degree and near-Boolean games are also finally discussed. 

Keywords: Set packing, Pseudo-Boolean function, Polynomial multilinear 
extension, Local search. 

MSC 2010: 06E30; 06A07; 06D72; 94C10. 

1 Introduction 

Consider a finite set N = {1,..., n} of items to be packed into feasible subsets, 
where these latter constitute a family T C 2^ = {A : A C N}. The problem 
is to find a subfamily J^* C T oi pairwise disjoint feasible subsets with largest 
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size IJ"*!. In the weighted version, a function w : T ^ K+ identifies as opti¬ 
mal those such subfamilies T* with maximum weight 

Maximizing |J^*| is equivalent to setting w{A) = 1 for all A ^ T. Thus this 
work proposes to use the polynomial multilinear extension, or MLE for short, of 
set functions (such as w) in order to evaluate families of fuzzy feasible subsets. 
Although unfeasible, such families shall still drive the search towards locally 
optimal feasible ones. 

Set packing is a key combinatorial optimization problem m extensively 
studied in computational complexity, where the aim is to find efficient algorithms 
whose output approximates optimal solutions within a provable bounded factor. 
In that field, the focus is placed mostly on non-approximability results for fc-set 
packing, where the size of every family member is no greater than some fc <g; n 
(and with unit weight for each member as above, [23]). Recall that if all family 
members have size k = 2, then the problem is to find a maximal matching 
in a graph with vertex set N, and an efficient (i.e. with polynomial running 
time) algorithm capable to output an exact solution is known to exist [I6j . In 
fact, if fc > 2, then k-set packing may be rephrased in terms of vertex clouring 
in hypergraphs, with special focus on the d-regular and /c-uniform case, where 
every element of the ground set is present in precisely d > 1 feasible subsets 
(i.e. |{A : i G A € J^}| = d for every i G N), each of which, in turn, has size k 
(i.e. |A| = k for every A G T) [Ilj . 

Set packing aslo has important applications, among which combinatorial 
auctions constitute a main and lucrative example: the ground set may consist of 
items to be sold in bundles (or subsets) towards revenue maximization, and once 
bids are processed the issue may be tackled as a maximum-weight set packing 
problem, with maximum received bids on bundles as weights [20) . Given the 
exponentially large size of the search space, revenue maximization often leads 
to use heuristics with no worst-case guarantee or, more simply, to sell each item 
independently but simultaneously over a sufficiently long time period El- 

The approach to set packing problems proposed in the sequel replaces stan¬ 
dard pseudo-Boolean optimization [3] with a novel near-Boolean method. While 
the former employs MLE to switch from {0,1} to [0,1] as the domain of each of 
the n variables, the proposed near-Boolean method relies on n variables ranging 
each in the 2"“^-set of extreme points of a unit simplex, and employs MLE to 
include the continuum provided by the whole simplex. The n variables corre¬ 
spond to the elements i G N oi the ground set, while the extreme points of each 
simplex are indexed by those subsets where each element is included. Then, 
the MLE of the resulting near-Boolean function allows to evaluate collections of 
fuzzy subsets of N or, equivalently, fuzzy subfamilies of feasible subsets A G J-. 
The objective function to be maximized takes thus values over the n-product of 
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2"’“^ — 1-dimensional unit simplices, allowing to design a flexible gradient-based 
local search. 

The following section comprenshively details the framework for the full¬ 
dimensional case T = 2^. This not only seems useful for applications, but most 
importantly allows to clearly see next that by simply introducing the empty 
set 0 and all n singletons {i} G 2^ into the family T of feasible subsets (with 
null weights ?u(0) = 0 = w({*} if {i} ^ J^) the whole class of set packing prob¬ 
lems may be handled by the proposed method. The gradient-based local search 
differs when switching from the full-dimensional case to the lower-dimensional 
one C 2^, in that with the latter a cost function c : —>• N also enters the 
picture, in line with greedy approaches to weighted set packing The cost 
c{A) = \{B : B G J", i? n A ^ 0}| of including a feasible subset in the pack¬ 
ing is the number of members with which it has non-empty intersection (itself 
included, hence c{A) G N, A G . 

Note that maximum-weight set packing may be tackled through constrained 
maximization of standard pseudo-Boolean function v : {0, —>■ K.+ 

given by u (a;^i,..., ) = X! xA;,w{Ak) , s.t. 

l<k<\J^\ 

Afc n A; 7 ^ 0 ^ XAk + XAi < 1 for all 1 < fc < Z < |J^|, 

where xa G {0,1} for aW A G F = {Ai,... A^jr^}. Also, v can be replaced 
with xMx ~ V, where M is a suitable |-F| x | J^|-matrix and x = {xai , ■ • ■, XA^J^^ ) 
[2] . An heuristic then finds a constrained maximizer x* , while the corresponding 
solution is F* = {A : x\ = 1}. This differs from what is proposed here, in many 
respects, the most evident of which being that v has |J^| constrained Boolean 
variables, while the expanded MLE developed below has n unconstrained near- 
Boolean variables. 


2 Full-dimensional case 

The 2"-set (0,1}" of vertices of the n-dimensional unit hypercube [0,1]" cor¬ 
responds one-to-one to the (power) set 2^ of subsets A C N through charac¬ 
teristic functions xa '■ N ^ {0,1}, A G 2^ defined by xa(*) = 1 if i G A and 
XAii) = 0 if f G N\A = A°, while collection {C(A, •) : A G 2^} is a linear basis 
of the vector space of real-valued functions w on 2^, where zeta function 
: 2^ X 2^ — >■ R is the element of the incidence algebra [mu] of Boolean lattice 
(2^, n, U) defined by C(A, B) = 1 ii B D A and C(A, B) = 0 if i? ^ A. Linear 
combination w{B) = X]yiG 2 " M'"(^)C(A, B) = B G 2^ applies 

to any w, with Mobius inversion fjF : 2^ —^ R uniquely given by (c is strict 
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inclusion) = 

= E (with CiB^A) = (-l)l^\^l) 

BCA 

= wlA) — E ^A’{B) (recursion, with ic(0) = 0). 

BcA 

Given this essential combinatorial “analog of the fundamental theorem of the 
calculus” [18], the MLE /*" : [0,1]" —^ R of w takes values w{B) = 


= rixB) 


Ag2" \ieA ) A<ZB 


on vertices, and f”’{q) = E (1) 

^ 62 ^^ VigA / 

on any point q= (gi,..., g„) G [0,1]". Conventionally, Hiefl •= 1 O P- 157]. 

Let = {A ■. i € A € 2^} = {Ai, ..., A 2 n-i} be the 2”“^-set of subsets 
containing each i G N. Simplex 


A, = 


^ 271 -] 

^)GRr': E 9f''=i] 

1 

1 


l<fe<2'*-i J 


has dimension 2" ^ — 1 and generic point qi € Ai. 

Definition 1 A fuzzy cover q specifies a membership distribution for each i G N 
over the 2"“^ subsets containing it, i.e. q= (qi, , qn) G A^v = Xi<i<„ A^. 

Equivalently, q = ^ A € 2^, q^ G [0,1]"} is a 2" — 1-set whose elements 

q'^ = {q^,..., q^) are n-vectors corresponding to non-empty subsets A G 2^ 
and specifying a membership for each i G N, with qf- G [0,1] if z G A while 
g^ = 0 if z G A“. Fuzzy covers being collections of points in [0, Ij", and the 
MLE /*“ of w allowing precisely to evaluate such points, the global worth lE(q) 
of q G Atv is the sum over all q^, A G 2^ of /“(g"^) as defined by (1). That is. 


hh(q) = E = E 






E 


BQA 

Vies / 


or lE(q) = E 

E 


Ae2^ 

BDA 

\i€A / 


p-{A). 


( 2 ) 


Example 2 For N = {1,2,3}, define w;({l}) = w;({2}) = ic({3}) = 0.2, 


4 






w({l,2}) = 0.8, w({l,3}) = 0.3, w({2,3}) = 0.6, w{N) = 0.7. Membership 
distributions of elements i = 1,2,3 over are gi € Ai, g 2 G A 2 , <73 € A 3 , 




' qk ' 


f \ 

^12 


„12 


„13 

Qi 


92 


93 


; 92 = 

qf 

, 93 = 

„23 

93 

V J 


\ J 


l 93^ / 


If q\^ = < 72 ^ = 1, then any membership (73 G A 3 yields 

W{qi,q2,q3) = w({l,2}) 

= ^i'({l,2}) + w({3}) = 1. 

This means that there is a eontinuum of fuzzy covers achieving maximum worth, 
i.e. 1. In order to seleet the one q = (<71,<72,<73) where = 1, attention must 
be placed only on exaet ones, defined hereafter. 

For any two fuzzy covers q = {q"^ : 0 ^ A € 2^} and q = {q^ : % ^ A ^ 2^}, 
define q to be a shrinking of q if there is a subset A, with > 0 and 


BCA 

In words, a shrinking reallocates the whole membership mass > 0 from 

A G 2^ to all proper subsets B C A, involving all and only those elements i € A 
with strictly positive membership q^ > 0 . 

Definition 3 Fuzzy eover q G A^v is exaet as long as W{q) ^ W{(i) for all 
shrinkings q of q. 

Proposition 4 If q is exaet, then | {i G A : < 7 ^ > O} | G {0, |A|} for all A G 2^. 
Proof: For 0 C ^^(q) = {i '■ qf > 0} C A, with a = |A+(q)| > 1, note that 

f^{q^) =J2 b<za+{ci) 

Let shrinking q, with q^ = q^ if B' ^ satisfy conditions 

1) X/B£2^n2'^’''(‘5) ~ di X/i3g2"n2'^+(‘5) ^ and 

2) n*eB = n*eB + n*eB for all B G 2 ^* A) such that |S| > 1. 


qf ABg A 


for all B €2^,ieN, 


[ 0 A B = A 

= qf+T. 9? fo" all z G ki. 

BcA 


5 



These are 2“ —1 equations with J2i<k<a ^ik) ^ variables qf, B C >l+(q), 
i G B. Thus there is a continuum of solutions, each providing precisely a 
shrinking q where 




This entails that q is not exact. ■ 

Partitions P = {Ai,... , A|p|} C 2-^ of are families of pairwise disjoint 
subsets called blocks [T], that is Ak Ci Ai = 0,1 < k < I < |P|, with union 
N = Ui<fc<|p| Ak- Any P corresponds to the collection {xa '■ A £ P} oi those 


|P| hypercube vertices identified by the characteristic functions of its blocks 
(see above). Partitions P can also be seen as p £ Ajv where pf = 1 for all 


A £ P,i £ A, i.e. exact fuzzy covers where each i £ N concentrates its whole 
membershisp on a unique A £ 2^, thus justifying the following. 

Definition 5 Fuzzy partitions are exact fuzzy covers. 

Ojective function W : A^v —t R includes among its extremizers (non-fuzzy) 
partitions. This expands a basic result in pseudo-Boolean optimization. Denote 
by ex{Ai) the 2"“^-set of extreme points of Ap For q £ Ajv,i £ N, let 
q = (?i|q_p with qi £ Ai and q_j e Ajv\i = ^jeN\i Aj. Then, for any w, 


H^(q) = E + E 

A'iz2N\rin 


Ae2f A'&2^\2t 



at all q £ An and for all i £ N. Now define 
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yielding W{q) = 


(3) 


Proposition 6 For all qG An, there are q,qG An 


such that 


(i) W{q) <W{q) < W{q) and, 


(a) q^, € ex{Ai) for all i G N. 

Proof: For i G N,q_^ G A^\j, define Wq_. : —>• IR by 


E n9/U“(i?uz). (4) 

BCA\i \jeB ) 

Let A+ = argmaxwq_^ and A“ = argminwq_., with A+ . 0 7 ^ A” at 

all q_j. Most importantly, 

Wi{q^\qi)= E ( 9 *^ (5) 

y4e2f 

where (•, •) denotes scalar product. Thus for given membership distributions of 
all j G N\i, global worth is affected by i’s membership distribution through a 
scalar product. In order to maximize (or minimize) W by suitably choosing qi 
for given q_j, the whole of z’s membership mass must be placed over A+ (or 
A“ .), anyhow. Hence there are precisely |A+ | > 0 (or |A“ . | > 0) available 
extreme points of . The following procedure selects (arbitrarily) one of them. 
RouNDUp(u;,q) 

Initialize: Set t = 0 and q(0) = q. 

Loop: While there is a z G with qi{t) ^ ex{Ai), 
set t = t + 1 and: 

(a) select some A* G A)^ 

(b) dehne, for all j G N, A G 2^ , 

f 9/(i-l)ifjV* 
qf{t) = \ 1 if j = z and A = A* . 

I 0 otherwise 


Output: Set q = q(t). 

Every change qf{t — 1) qf{t) = 1 (for any z G N,A G 2f) induces a 
non-decreasing variation W{q{f)) — W{q{t — 1)) > 0. Hence, the sought q is 


7 



provided in at most n iterations. Analogously, replacing A+ with yields 
the sought minimizer q (see also [3l p. 163]). ■ 

Remark 7 For i G N, A G , if all j G A\j ^ 0 satisfy = 1, then (4) 
yields Wq_.{A) = w(A) — w{A\i), while u;q_^({j}) = r/;({j}) regardless of 

Corollary 8 Some partition P satisfies W{p) > W{q) for all q G Ajv, with 


Wip) = Eagp^(^)- 


Proof: Follows from propositions 4 and 6 (with the above notation associating 
p G Ajv to partition F). m 

Defining global maximizers is clearly immediate. 

Definition 9 Fuzzy partition q G Ajv is a global maximizer ifW{q) > W{q) 
for all q G A^. 

Concerning local maximizers, consider a vector w = (a;i,...,a;„) G R++ of 
strictly positive weights, with wjv = focus on the (Nash) equi¬ 

librium [13] of the game with elements i G N a.s players, each strategically 
choosing its membership distribution qi G A^ while being rewarded with frac¬ 
tion ■^W{qi, ..., q„) of the global worth attained at any (qi,..., q„) = q G Ajy. 

Definition 10 Fuzzy partition q G Ajv is a local maximizer if for all qi G A^ 
and all i G N inequality Wi{qi\q_i) > Wi{qi\q_i) holds (see (3)). 

This definition of local maximizer entails that the neighborhood A/'(q) C A^v of 
any q G A^r is 



Definition 11 The {i, A)-derivative of W at qG A^r is 


dW{q)/dqt = W{q{l,A)) - W{q{l,A)) 





q^ for all j G N\i, B G 2^ 


1 for j = i, B = A 
0 for j = i,B A 
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andq{i,A)= A),... A)j given by 

/ '?f V all j e N\i, B € 2f 
’ ^0 for j = i and all B G 

thus yW{q) = {dW{q)/dqf : i € N, A G 2^} G K”^” ^ is the (full) gradient 
ofW at q. The i-gradient Viiy(q') € ofW at q— qi\q_i is set function 

^iW{q) : > R defined by '\/iW{q){A) = Wq_.{A) for all A G 2f, where 

Wq_. is given by ()). 

Remark 12 Membership distribution qf^i,A) is the null one: its 2"“^ entries 
are all 0, hence qff,A) ^ A^. 

The setting obtained thus far allows to conceive searching for a local maximizer 
partition p* from given fuzzy partition q as initial candidate solution, and while 
maintaing the whole search within the continuum of fuzzy partitions. This idea 
may be specified in alternative ways yielding different local search methods. 
One possibility is the following. 

LocalSearch(w, q) 

Initialize: Set t = 0 and q(0) = q, with requirement |{z : q(^ > 0}| G {0, |Gl|} 
for all A G 2^. 

Loop 1: While 0 < Yl,ieA ^fi^) < 1^1 for a A G 2^, set t = t + 1 and 

(a) select a A*{t) G 2^ such that 


ieA*{t) jeB 

for all B G 2^ such that 0 < < l-®l> 

(b) for i G A*{t) and A G 2f , define 


qfit) 


1 AA = A*{t), 
Oil A j^A*{t), 


(c) for j G N\A*{t) and A G 2^ with Ar\A*{t) = 0, define qf{t) = qf{t — l)+ 



/ \ 


1 \ 

+ 

wiA) 


E 


BG2f 


B'e2f 


y BnA*(t)7i0 y 


yB'nA*(t)=0 j 
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(d) for j G N\A*{t) and A G 2^ with A n A*{t) ^ 0, define 


g/(t)=0. 


Loop 2: While qf{t) = 1, |A| > 1 for a i € and w{A) < w({i}) + w{A\i), set 
t = t + 1 and define: 




1 if \A\ = 1 

0 otherwise 

1 if S = A\i 
0 otherwise 


for all i G 2f, 


for all j G B G 2 


N 

j ’ 


1 ) for all j' GA‘=,BG2f,. 


Output: Set q* = q(t). 

Both RoundUp and LocalSearch yield a sequence q(0),..., q(t*) = q* 
where q* G ex(Ai) for all i G N. In the former at the end of each iteration t the 
novel q(t) G Af{q{t — 1)) is in the neighborhood of its predecessor. In the latter 
q(t) ^ Af{q{t — 1)) in general, as in |P| < n iterations of Loop 1 a partition 
{^*(1),..., A*(|P|)} = P is generated. Selected subsets or blocks A*{t) G 2^, 
t = 1,...,|P| are any of those where the sum over members i G A*{t) of 
(i, A*(t))-derivatives 9W(q(t —l))/9(;f ^*\t—l) is maximal. Once a block A* (t) 
is selected, then lines (c) and (d) make all elements j G N\A*{t) redistribute 
the entire membership mass currently placed on subsets A' G 2^ with non¬ 
empty intersection A' n A*{t) ^ 0 over those remaining A G 2^ such that, 
conversely, A n A* it) = 0. The redistribution is such that each of these latter 
gets a fraction w{A)/^g^ 2 ^-BnA*{t)=<li''^iB) of the newly freed membership 
mass X]A'e 2 "-A'nA*(t) 5^0 ~ !)• subsequent Loop 2 checks whether the 

partition generated by Loop 1 may be improved by exctracting some elements 
from existing blocks and putting them in singleton blocks of the final outpnt. 
In the limit, set function w may be such that for some element i G N global 
worth decreases when the element joins any subset A G 2f^, |A| > 1, that is to 
say w(A) - wiA\i) - u;({i}) = 'Eb(^ 2 ^\ 2 ^\^:\b\>i < 0- 

Proposition 13 LocalSearch(IT, q) outputs a local maximizer q*. 

Proof: It is plain that the output is a partition P or, with the notation of 
corollary 8 above, q* = p. Accordingly, any element i G N is either in a 
singleton block {i} G P or else in a block A G P, i G A such that |A| > I. 

|i| 

In the former case, any membership reallocation deviating from p] ‘ = 1, given 
memberships pj,j G N\i, yields a cover (fuzzy or not) where global worth is the 
same as at p, because Y\j(sB\iPf ~ ^ ^ 2 (^\A (see example 2 above). 
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In the latter case, any membership reallocation qi deviating from pf = 1 (given 
memberhips pj,j € N\i) yields a cover which is best seen by distinguishing 
between 2^\A and A. Also recall that w{A) — w{A\i) = 

Again, all membership mass ^Be 2 ^\A if > ^ simply collapses on singleton {i} 
because = 0 for all B G 2f^\A. Hence W{p) — W{qi\p_^) = 


= w{A) - u;({i}) - qf 




p^(B) 


E 


pr{B') = 


= {pf - <lf) 


Be2'^\2"i\G|B|>l 


Now assume that q is not a local maximizer, i.e. lH(p) — W{qi\p_j) < 0. Since 
pf — qf > 0 (because pf = 1 and qi G is a deviation from pi), then 

pA{B) = w{A) — w{A\i) — icdi}) < 0. 

Hence q cannot be the output of Second Loop. ■ 

In local search methods, the chosen initial canditate solution determines 
what neighborhoods shall be visited. The range of the objective function in a 
neighborhood is a set of real values. In a neighborhood Af{p) of a p G Ajv or par¬ 
tition P only those X]y 4 eP-|A|>i 1^1 elements i G A in non-sigleton blocks A £ P, 
|A| > 1 can modify global worth by reallocating their membership. In view of 
(the proof of) proposition 13, the only admissible variations obtain by deviat¬ 
ing from pf = I with an alternative membership distribution qf G [0,1), with 
W{qi\p_i) -W{p) equal to {qf - I) EBe 2 ^\ 2 Mi,|B|>i P^{B) + (1 - qf )w{{i}). 
Hence, choosing partitions as initial candidate solutions of LocalSearch is 
evidently poor. A sensible choice should conversely allow the search to ex¬ 
plore different neighborhoods where the objective function may range widely. 
A simplest example of such an initial candidate solution is qf = for all 

A £ and all i £ N, i.e. the uniform distribution. On the other hand, the 
input of local search algorithms is commonly desired to be close to a global 
optimum, i.e. a maximizer in the present setting. This translates here into the 
idea of defining the input by means of set function w. In this view, consider 
qf = w(A)/X;Bg 2 f ^(B), yielding ^ ^ for all A,B £ 2^0 2^ and 

all i, j £ N (see lines (c), (d)). 

With a suitable initial candidate solution, the search may be restricted to ex¬ 
plore only a maximum number of fuzzy partitions, thereby containing the com¬ 
putational burden. In particular, if q(0) is the finest partition {{!},..., {n}} 
or ql''\o) = 1 for all i £ N, then the search does not explore any neighbor- 
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hood at all, and such an input coincides with the output. More reasonably, let 
A^ax = {Ai ,..., Afe} denote the collection of D-maximal subsets where input 
memberships are strictly positive. That is, qf'^' > 0 for all i G 1 < fc' < fc as 
well as < 7 ^ = 0 for all B G 2^\ (2^^ U • • • U and all j G B. Then, the output 
shall be a partition P each of whose blocks A G P satisfies A C Ak' for some 
1 < k' < k. Hence, by suitably choosing the input q, LocalSearch outputs a 
partition with no less than some maximum desired number fc(q) blocks. 

3 Lower-dimensional case 

If -F C 2^^, then 2f^ n -F 7 ^ 0 for every i G N, otherwise the problem reduces to 
packing the proper subset N\{i : A" fl 2f^ = 0} of elements contained in at least 
one feasible subset. As outlined in section 1 , without additional notation simply 
let {0} G T 5 {d for all i G with null weigths w(0) = 0 = w({i}) if {d ^ 
Thus {P, A) is a poset (partially ordered set) with bottom element 0, and poset 
function w : P ^ R+ has its Mobius inversion fp : P ^ M. [18]. Memberships 
qi distribute over P^ = 2f r\ P = {Ai,..., A|_ 7 r.|}, with lower(|Ai|)-dimensional 
unit simplices 



and corresponding fuzzy covers q G = Xi<i<ji A^. Note that a fuzzy cover 
now may maximally consist of |.F| —1 points in the unit n-dimensional hypercube 
[0,1]". Hypercube [0,1]" is replaced with C{P) = co{{xa ■ A G J^}) C [0,1]”, 
i.e. the convex hull of feasible characteristic functions, regarded as n-vectors 
[2. Recursively (with w(0) = 0), Mobius inversion ■. P is 

^P{A)=w{A)- Y. 

bgt-.bga 

while the MLE /“ : C{P) —>• K. of ui is 

r{q^)= E 

BeBn2'^ \i&B / 

Therefore, every fuzzy cover q G A^v has global worth 

^(q) = E E 

AgB BGBn2^ \ieB ) 
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For all i £ N,qi G A*, and q_i € A^\ j = x Aj 

jGN\i 


W_.(q^,) = ^ 

E 

\ n cn 

AGBi 


) 


+ E 

E 1 

(n 9/1 GiB') 


_B'6JPn2^' 

j J 


W^iq^\(l^i) = E <lf 

E 

i n if] 

AGBi 

BeBin2'^ 

\jeB\i j 


yielding again 

W^(q) = W^*fe|q-J + VF_,(q_J. (6) 

From (4) above, Wq_. : —5- R now is 


Wq^M)= E n (7) 

-B6^in2'^ \jeB\i J 

for all i G N, all A G Ti and all q_^ € A^\j. 

For each i G N, denote by ex(Aj) the set of \Fi\ extreme points of simplex 
Ai. Like in the full-dimensional case, at any fuzzy cover q € Ajv every i G N 
such that Qi ^ ex{Ai) may deviate by concentrating its whole membership on 
some A G Ti such that r(;q_.(A) > Wq_.{B) for all B G Ti. This yields a 
non-decreasing variation fF((j'i|q, -*) > LF(q) in global worth, with qi G ex{Ai). 
When all n elements do so, one after the other while updating rCq_.(t) as in 
RoundUp above, i.e. t = 0,1,..., then eventually q = {qi,...,qn) is such 
that q G X ex{Ai). Yet cases T C 2^ and F = 2^ are different in terms of 

ieN 

exactness. Specifically, consider any % ^ A G F with |{i : qf = 1}| ^ {0, |A|} 
or GI+ = {* : qf = 1} C A, with 

is likely to admit no shrinking (see above) yielding an exact fuzzy cover with 
same global worth as (non-exact) q. 

Proposition 14 The values taken on exact fuzzy covers do not saturate the 
range of W : An -g 1R+. 

Proof: By example: N = {1,2, 3,4} and F = {A, {4}, {1, 2}, {1,3}, {2, 3}}, 
with worth w{N) = 3, ri;({4}) = 2, t(;({i, j}) = 1 for 1 < i < j < 3. Now 
define q = (gi,..., < 74 ) by = 1 = q^, i = 1,2, 3, with non-exactness due to 
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inequality |{j : > 0}| = 3 < 4 = |A^|. As 


M^(q) = u;({4}) + ^ w;({i, j}) = 2 + 1 + 1 + 1 

l<i<j<3 


and A+ = {1,2,3}, for all distributions < 71 ,( 72,93 placing membership only over 
feasible B £ F f] global worth is 11 ^( 91 ,( 72 , 93 ,( 74 ) < ll^(q). ■ 

This simple result is useful because it indicates that an inaccurate search 
among optimal fuzzy covers may lead to a maximizer, either global or local, 
which is not reducible to any feasible solution of the original set packing prob¬ 
lem. In the present setting, such feasible solutions are partitions P all of whose 
singleton blocks with worth 0 are not included in the packing. In fact, simi¬ 
larly to the full-dimensional case, fairly simple conditions may be shown to be 
sufficient for a partition to be a local maximizer. 

Definition 15 Any qi\q_i = q & An is a local maximizer of W : An ^ R+ if 
Wi{qi\q-i) A kki(9i|9-i) for all i G N and all qi G Ai (see (6) above). 

The neighborhood A/’(q) C A^r of q G Ajv thus is 

■^(q) = U {q ■ q = ^dq-oi* s a,|. 

ieN 

Evidently, there are many partitions P with associated p such that p G An 
( with the notation of corollary 8 above). For example P± = {{I},...,{n}}, 
i.e. the bottom element of the geometric lattice (7^'^,A,V) of partitions of JV 
ordered by coarsening [T] . Other simple examples are given, for every A G P, hy 
the corresponding modular element [32] {A} U of (P^, A, V) whose unique 
non-singleton block is A. In fact, any partition P such that A G P for each 
block A G P has associated p satisfying p G A^r. 

Proposition 16 Any partition P with associated p such that p G An is a local 
maximizer if for all A G P 

w{A) > w({f}) -P iP’iB). 


Proof: Firstly note that for all blocks A G P, if any, such that |A| = 1 there 
is nothing to prove, as the summation reduces to w(0) = 0, and thus there 
only remains w({f}) > w({f}). Accordingly, let A G P and |A| > 1. For every 
i G A, any membership reallocation qi G A^ deviating from pi (i.e. p'^ = 1), 
given memberships p_j of other elements j G N\i (i.e. A^ ^ p^ = 1 for all 
A' G P and all j G A'), yields q = (9i|p_j) G An which is best analyzed by 
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distinguishing between Fi\A and A. In particular, 


n;(A) =«;({*})+ ^ ^ 

| B |>1 


All membership mass > 0 collapses on singleton {i}, because 

Y\ pf> = ^ for all B e Ti\A by the definition of p_j (see example 2 above). 

i' ^B\i 

Thus, 


lT(p) -W{q^\Tp_i) = w{A) -icdiDT 

( \ 

qf E 

BGBin2^ BGBn2^\^ 

\ \B\>1 j 

= {pf-Qf) E 

Be.Fin2'^ 

|B |>1 


Now assume that p is not a local maximizer, i.e. IT(p) — Il^(gi|p_j) < 0. Since 
pf — qf > 0 because pf = 1 and qt € is a deviation from pi , then 

^ ^»(B)=^A)-n;({*})- ^ p^{B)<0 

|B |>1 


must hold. This contradicts the premise w{A) > w({i}) + BBGJ^n 2 ^\* 
for all A G P and i G A, thus completing the proof. ■ 

4 Local search with cost 

In order to design a gradient-based local search for this lower-dimensional case, 
the only tool still missing is the derivative, which clearly shall reproduce defi¬ 
nition 11 above with Pi in place of 2^. Before that, as outlined in section 1, 
let c : P N count the number c{A) = \{B : B G B, B D A ^ 9}\ of feasible 
subsets with which each A G B has non-empty intersection, itself included, i.e. 
c{A) G {1,..., Ip’ll is the cost of including A in the packing. Accordingly, the 
underlying poset function w : B ^ M+ now used (still taking positive values 
only) incorporates both weights w{A),A G B (used thus far) and costs by means 
of ratio w{A) = 7 ^- The result is quasi-objective function W : An R+, 
obtained via MLE /“ : C{B) -G IR+ of w, i.e. IE(q) = 
of the above applies invariately simply replacing W with W. 
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Definition 17 The {i, A)-derivative of W at q€ Ajv, 


A e J-., IS dW{q)/dqt = W{q{l,A)) - W{q{l,A)) = 


= W^{qS, A)\q_S, -W^[qf^i,A)\q_.{i,A)^, 

withq{i,A)= A),... A)j given by 

( q^ for all j G N\i, B G Tj 
lforj=i,B = A 
[ 0 for j = i, B ^ A 


andq{i,A)= (^q^{i, A),... ,q^{i, A)j given by 




q^ for all j G N\i, B G Tj 
0 for j = i and all B G Ti 


The (full) gradient of W at q G is 

yw{q) = ^dW{q)/dqf : i G N, A G G l-^^l 

as well as the i-gradient ViW{q) G of W at q = {qi\q^i) G Ajq is poset 

function ViW{q) : Bi ^ M. defined by 'ViW{q){A) = Wq_.{A) for all A G Bi, 
where Wq_. is given by (7) with w in place of w. Again, membership distribution 
q^{i,A) is the null one: its \Bi\ entries are all 0, hence q^{i,A) ^ A^. 


LS-WithCost(w, q) 

Initialize: Set t = 0 and q(0) = q, with requirement \{i : qf > 0}|G{0,|GI|} 
for all A G B, w{A) > 0. 

Loop 1: While 0 < ^ti^) < 1^1 for a A G set t = t + 1 and: 

(a) select a A*{f) G B such that 


i^A*{i) j€B 

for all B G 2^ such that 0 < ^ 1^1’ 

(b) for i G A*(t) and A G Bi, define 




1 iIA = A*{t), 
Oil A ^A*{t), 
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(c) for j G N\A*{t) and A € Tj with Af\A*{t) = 0, define qf{t) = qf{t — \)+ 


\ 


w{A) Y. 


V 


BnA*{t)=^(b 


) 




E 


\B' r\A* {t) = 


) 


(d) for j G N\A*{t) and A G J^j with A n A*{t) ^ 0, define 


qtit)=0. 


(e) for A G with An A* (t) =0, update cost function by 

c{A) = \{B : B eT,BnA^9 = Bn A*(t)}| 
and plug it into w. 

Loop 2: While qf(t) = 1, |A| > 1 for a i G fV and 
w{A) < wdi}) + E 


set t = t + 1 and define: 


qfit) 

qf'{t) 


1 if |i| = 1 
0 otherwise 

1 if B = A\i 
0 otherwise 


for all A G Bi, 
for all j G A\i, B G Bj, 


qf,{t - 1) for all j' gA^,Bg By. 


Output: Set q* = q(t). 

Both LocalSearch and LS-WithCost generate in |P| < n iterations of 
Loop 1 a partition {A*(l),..., A*(|P|)} = P. Now selected blocks A*{t) G B, 
1 £ ^ < |-P| are any of those feasible subsets where the minimum over elements 
i G A*{t) of (i, A*(t))-derivatives dW{c\{t — \))ldqf — 1) is maximal. The 
following Loop 2 again checks whether the partition generated by Loop 1 may 
be improved by exctracting some elements from existing blocks and putting 
them in singleton blocks of the final output, which thus allows for the following. 


Proposition 18 LS-WithCost(W, q) outputs a local maximizer q*. 

Proof: Follows from proposition 16 since Loop 2 deals with w, not with w. ■ 
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Concerning input q = q(0), consider again setting 


<7^ = 


)(A) 


w{B) 


for all A G -Fi, i € iV, which entails 


4 = = % i„ all A,Be € «■ 

q° w{B)c[A) qf 

Evidently, Loop 1 may take exactly the same form as in LocalSearch, 
that is with selected blocks A*{t) G T,t = 1, • ■ •, 1^*1 of the generated partition 
P being any of those feasible subsets where the sum, rather than the minimum, 
over elements i G A*{t) of (i, ^*(t))-derivatives dW{c{{t — l))/dqf — 1) is 
maximal. This possibility seems appropriate in applicative scenarios, where set 
packing is mostly dealt with in its weighted version. Yet, using the minimum 
in place of the sum, although computationally more demanding, appears in¬ 
teresting for fc-uniform set packing problems (see section 1), widely studied in 
computational complexity. In fact, for the A:-uniform case Mobius inversion is 
pP’{A) = if |yl| = k and fjP{A) = 0 if |A| G {0,1} for &\\ A G T (recall the 
convention {0} G T ^ {i} for all i G N), with the cost function updated at each 
iteration according to line (e). It is also evident that in A:-uniform set packing 
Loop 2 is ineffective. 


5 Near-Boolean functions 

Boolean functions [5] provide key analytical tools and methods with a variety 
of important applications. Beyond set packing problems that here constitute 
the main benchmark, this section further develops the full-dimensional case de¬ 
tailed in section 2 with the aim to indicate additional opportunities obtained 
from expanding the standard framework where pseudo-Boolean models are tra¬ 
ditionally exploited. Recall that Boolean functions of n variables have form 
/ : {O,!}" —>■ {0,1}, and constitute a subclass of pseudo-Boolean functions 
/ : {0,1}” —>■ R, which in turn admit the unique MLE / : [0,1]" —>■ R over the 
whole n-dimensional unit hypercube extensively employed thus far. The n vari¬ 
ables thus range each in the unit interval [0,1]. Such a setting is here expanded 
by letting each variable i = 1 ,... ,n range in a 2 "“^ — 1 -dimensional simplex 
Ai , with the goal to evaluate collections of fuzzy subsets of a n-set through the 
MLE given by (1) and (2). 

Definition 19 Near-Boolean functions of n variables have form 

F : X ex{Ai) —> R. 

l<2<n 
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Following [13 p. 4], denote by N = {l,...,n} the set of indices of variables 
(i.e., the ground set in previous sections). As already observed, any pseudo- 
Boolean function has a unique expression as a multilinear polynomial / in n 
variables: f{xi,... ,x„) = ^ I a a 0 ) ■ Uniqueness is customarily shown 

by induction on the size 0 < |A| < n of subsets of variables 13 p. 162], although 
it seems a by-product of the Mobius inversion of (po)set functions [18], which is 
unique indeed, i.e. A G simply is the Mobius inversion of some unique 
set function w ■. 2^ such that w{A) = /{xa), where XA is the characteristic 
function defined in section 2, i.e. y : 2^ —>• {0,1}" with x(^) = XA- 

Definition 20 The MLE F of near-Boolean functions F has polynomial form 


F : X Ai 

l<i<n 


given by expression (2) in section 2, that is 


Hq) = E 


E 


BDA 

\ieA / 




with (see above) q— (gi,..., g„) and qi = {qf^ ,..., qf'^"' ^€ A. 


5.1 k-degree approximations 

In line with uni, the issue of approximating a given near-Boolean function F 
by means of the least squares criterion amounts to determine a near-Boolean 
function Fk such that 


E [F(q)-F,(q)]^ (8) 

qG X ex{Ai) 


attains its minimum over all near-Boolean functions Fk with polynomial MLE 
Fk of degree k, that is 


A(q) = E 


!Aj<fe 


E 


BDA 

\iGA / 




or, equivalently stated in terms of the underlying set function ic, such that 
^“(A) = 0 if |A| > k. 

Near-Boolean functions F take their values on n-product x ex(Ai), and 

iGN 

|ex(Ai)| = 2"“^ for each i G N. They might thus be regarded as points 
F G in a 2"U-i)_£tiinensional vector space. In view of proposition 
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4 concerning exactness, this seems conceptually incorrect and with useless enu- 
merative demand. Specifically, for every partition P € with associated 
p S X ex(Ai), there clearly exist many non-exact q S x ex(Ai) such that 

ieN iGN 

F(q) = F{p) (see corollary 8 above). Counting these redundant extreme points 
of simplices appears wothless. For this reason, fc-degree approximation is dealt 
with by replacing expression (8) with the following, applying to partitions p 
only 

^ [F(p)-Ffc(p)f. (9) 

pG X ex{Ai) 

i€N 

The number \V^\ of partitions of a n-set is given by Bell number Bn [171 [1]. 
Accordingly, near-Boolean functions might be regarded as points F € K®" in 
a S„-dimensional vector space. Still, this also is far too large, as points in 
such a vector space correspond in fact to generic partition functions, i.e. with 
Mobius inversion free to live on every partition P G . Conversely, near- 
Boolean functions factually involve only partition functions h : -G R such 

that h{P) = h^{P) — for some set function ic : 2^ —>• R. The 

Mobius inversion of these partition functions lives only on the 2" — u modular 
elements [22] of lattice (7^'^,A, V), namely on those partitions with a number 
of non-sigleton blocks < 1. In turn, this entails that when regarded as points 
in a vector space (i.e. expressed as a linear combination of a basis, see above) 
these functions may be seen as h^, G R^"“”. This can be shown via recursion 
through the Mobius inversion of additively separable partition functions; it is 
here omitted being contained in mis]. 

When all these facts are properly taken into account, the issue of fc-degree 
approximation for near-Boolean functions is seen to reduce to the same issue 
for traditional pseudo-Boolean functions, which is already exhaustively detailed 
in nuj. What is crucial emphasizing though, is that while for pseudo-Boolean 
functions there exists a unique best fc-degree approximation for all 0 < fc < n, on 
the other hand every near-Boolean function admits a continuum of set functions 
w determining their unique best fc-degree approximation, and this applies to all 
0 < k < n. Furthermore, case fc = 0 is of no concern in that the emptyset cannot 
be a block of any partition. In particular, consider the linear case, which is the 
main case in [lOl p. 4]. The issue is to find a best (least squares) approximation 
Fi of any given F. That is, the set function w determining Fi has to satisfy 
u>{A) = ^ ^ Then, 

hw{p) = = 

AeP AePieA 

for all P G . Thus, hw is a constant partition function, i.e. a valuation 
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[T] of partition lattice (7^^,A,V). Also, any further linear v : 2^ ^ M. such 
that v{N) = w{N) satisfies hy{P) = hy,{P) for all P G . This means that 
there is a continuum (i.e. a non-unit n — 1 -dimensional simplex) of equivalent 
linear v ^ w such that hy, = hy, obtained each by distributing arbitrarily the 
whole of w{N) over the n singletons {i} G 2^. Cases A: > 1 are still all the 
same. To see this, consider a set function w such that yP{A) ^ 0 for one or 
more (possibly all (^)) subsets A g2^ such that |A| = k. Now fix arbitrarily n 
values u({z}), i G N with For all A G 2^, \A\ > 1 

Mobius inversion fp : 2^ ^ M. can always be determined uniquely through 
recursion by 

v{A) + ^ u({*}) = ^ ^ u({*}) = 

ieA” BCA ieA” 

= w{A) + w({i}) = Y ^({0)- 

B<ZA i^A'= 

One thing must be absolutely clear: there is a continuum of equivalent set 
functions (i.e. w and v) available for the sought k-degree approximation Fk, 
but still the Bn values taken by are unique and indpendent from the chosen 
set function in the continuum available (and thus such Bn values are also unique 
for any fixed F to be approximated). Any F clearly is its own best n-degree 
approximation. 

Remark 21 There is a continuum of set functions w equivalently determining 
the MLE F of F. 

5.2 Near-Boolean games 

In view of the above definition of local maximizers relying on equilibrium con¬ 
ditions for strategic n-player games, and having mentioned additive separablity 
of partition functions or global games [7], it seems now natural to consider 
variables as players in near-Boolean games (see also [B section 3]). 

Definition 22 A near-Boolean n-player game is a triple {N, F, tt) such that 
N = {1,... ,n} is the player set and F is a near-Boolean function taking real 
values on profiles q = (gi,...,g„) G x ex{Ai) of strategies, while payoffs 

ieN 

TT : X ex{Ai) -A M" are efficient, namely 7r(q) = ... ,'Kn{q)) satisfies 

ieN 

Ml) = F{q) at all q G x ex{Ai). 

ieN 

Definition 23 A fuzzy near-Boolean n-player game is a triple {N, F, n) such 
that N = {!,..., n} is the player set and F is the MLE of a near-Boolean 
function taking real values on strategy profiles q = (gi,..., g„) G x while 

ieN 
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TT : X Ai ^ R” efficiently assigns payoffs 7r(g) = (7ri(q),... ,TTn{Q)) to players, 
i£N 

ali x A^. 

ieN 

Game theory since many years is mostly concerned with finite sets of players, 
which is indeed the case both for near-Boolean games and fuzzy ones. Given 
this, a main distinction is between games where players have either finite or 
else infinite sets of strategies, with near-Boolean games in the former class and 
fuzzy ones in the latter. In addition, players may play either deterministic 
(i.e. pure) or else randomized (i.e. mixed) strategies. In the latter case equi¬ 
librium conditions are stated in terms of expected payoffs, and by means of 
fixed point arguments for upper hemicontinuous correspondences such condi¬ 
tions are commonly fulfilled [131 P- 260]. The sets of deterministic strategies 
in fuzzy near-Boolean games are precisely the sets of randomized strategies in 
near-Boolean games. Nevertheless, the payoffs for the fuzzy setting evidently 
are not expectations. 

The main framework where (possibly fuzzy) near-Boolean games seem ap¬ 
propriate is coalition formation, which combines both strategic and cooperative 
games. A generic strategy profile q G x ex{Ai) of near-Boolean (non-fuzzy) 

ieN 

games may well fail to be exact (see proposition 4), but it is understood at this 
point that there is a unique partition P of N with associated p G x ex{Ai) such 

ieN 

that i^(p) = T’(q). Let p(q) be such a unique p. In view of the above discussion 
on approximations, it is also clear that for every p there are many q such that 
p=p(q). In these terms, near-Boolean games model stategic coalition formation 
in a very handy manner, in that they totally by-pass the need to define a mecha¬ 
nism mapping strategy profiles into partitions (or coalition structures) of players 
[21] . More precisely, a mechanism is a mapping M : x ex(Ai) ^ such that 

ieN 

when each player i € N specifies a coalition Ai G 2)^, then M{Ai, ..., A„) = P 
is a resulting partition. If the n specified coalitions Ai,i G N are such that 
for some partition P it holds Ai = A for all i G A and all A € P, then 
M{Ai,..., An) = P. Otherwise, the partition P' = M(Ai,...,A„) gener¬ 
ated by the mechanism shall be a rather fine one, i.e. possibly consisting of 
many small blocks (depending on the chosen mechanism). Gonversely, near- 
Boolean games do not need any mechanism, in that even if players’ strategies 
(gi,..., Qn) = q are such that q does not correspond to a partition, still the 
global worth F{q) is that attained at the partition P with corresponding p(q), 
i.e. whose blocks A £ P each include maximal subsets of players choosing the 
same A' A A. 

Now fix a set function or coalitional game v : 2^ ^ R+, u(0) = 0 such that 
^(p) = partitions P, and let the payoffs be defined, for all 
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i ^ A and all A G P, by 




\B\ ’ 

BG2^\2'^\' 


where P is the partition with associated p(q). This is in fact a well-known 
coalition formation game, where payoffs are given by the Shapley value [IS]- 

Definition 24 A local maximizer q G x ex{Ai) of near-Boolean function F 

ieN 

satisfies for all i € N and all q[ G ex{Ai) inequality F{q) > F(g'|q_j). 

Remark 25 If payoffs are given by ni{q) = for all i G N, with 

wi ,... ,ujn >0, then near Boolean games are (pure) common interest potential 
games The set of equilibria of near-Boolean game (N, F, tt) coincides 

with the set of local maximizers of F. 


6 Conclusions 

Via polynomial MLE, the proposed near-Boolean functions of n variables take 
values on the n-product of high-dimensional unit simplices. This enables to 
approach discrete optimization problems, namely set packing, with an objective 
function defined over a continuous domain, with feasible solutions found at 
extreme points of the simplices. Approximations with polynomials of bounded 
degree are finally discussed, while near-Boolean n-player games provide a new 
modeling of strategic coalition formation. 
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